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New physics models that increase the decay rate of B^ — >■ t^t~ contribute to the absorptive 
part of Bs-Bs mixing, and may enhance AFs all the way up to its current experimental bound. In 
particular, the model with a scalar leptoquark can lead to a significant violation of the expectation 
AFs < AFs (SM). It can even allow regions in the ATs~Ps parameter space that are close to the best 
fit obtained by CDF and D0 through Bs — >■ J/tpcj>. In addition, it can help explain the anomalous 
like-sign dimuon charge asymmetry observed recently by D0. A measurement of BR{Bs — ^ t'^t~) 
is thus crucial for a better understanding of new physics involved in Bs-Bs mixing. 
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I. INTRODUCTION 

In the standard model (SM), the Cabbibo-Kobayashi- 
Maskawa (CKM) mixing matrix is the only source of 
charge-parity (CP) violation. The data from the decays 
of K, D and B mesons have so far been consistent with 
this paradigm, however the flavor changing neutral cur- 
rent (FCNC) processes involving 6 — >■ s transitions are 
expected to be sensitive to many sources of new physics 
(NP) [l|. This is why the B^ meson is one of the most 
important and interesting portals for indirect detection 
of such NP models. 

In this paper we shall concentrate on the oscillation 
parameters in the Bg-Bg system. The average de- 
cay width Ts = (TsH + Tsl)/2 = (0.6791°:°}?) PS"' 
and the mass difference AMg = Msh — Msl ~ 
(17.77 ± 0.10 ± 0.07) ps^^ have already been measured to 
an accuracy of better than '^ 2% [2-|3| and play an impor- 
tant role in constraining any new physics. Here the labels 
L and H stand respectively for the light and heavy mass 
eigenstates in the neutral Bs system. The decay width 
difference AFs = Tsl — Tsh and the Bg-Bg mixing phase 
are relatively less certain. The SM predictions for these 
quantities are [5| 

^^^ - (0.096 ±0.039) ps~^ 

VcbV*g\ -; 0.019 ±0.001, (2) 



AF^ 
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^J/^HSM) ^ Arg 
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where 2/3, 



J/i"l> 



is the mixing phase relevant for Bg 




decay. The recent CDF and D0 measurements 
using the angular analysis in Bg — )■ J/ip4' decay 
givelll 

(3) 
(4) 



AF, ^ ±(0.154l°:°fo) ps-i , 

pj,^, = (o.39t2:l«)u(i.i8l°:i^) 



where the second set in the last line is just the comple- 
ment of 7r/2 for the first set. This reflects the ambiguity 
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FIG. 1: The combined experimental constraints by CDF and 
D0 through Bs — >■ J /^<l>- Blue, red and green contours (from 
inner to outer) correspond to the 68%, 95% and 99% C.L. 
regions. The sinusoidal green band corresponds to the relation 
AFs ~ AFs cosi^s, valid when NP does not contribute to 

Fl2s. 



The figure is taken from Q. 



in the determination of (3s' . Note that the sign of AF^ 
is undetermined. The positive and negative signs corre- 
spond, respectively, to the two disconnected regions in 
the allowed parameter space for (3s . Alternative ways 
of removing this sign ambiguity have been suggested in 
[FLI ] . The correlated constraints are shown in Fig. [TJ The 
SM prediction for (AFs,/3s' ) is excluded by the data 
to 90% C.L.. Hence the exploration of new physics effects 
on these quantities becomes imperative. 

While many new physics models can affect (3's and 
make its value anywhere in its conventional allowed range 
[— 7r/2,7r/2], the ability of new physics to influence AFs 
is rather limited. Indeed, the width difference is 



AFs = 2|Fi2Jcos( 



(5) 



where (ps = Arg(-Mi2s/ri2s)- Here Mi2s and Ti2s are 
the dispersive and absorptive parts, respectively, of the 
Bg-Bs mixing amphtude. In the SM Q 

(/)« = 0.0041 ± 0.0007 , (6) 

and hence Arf'^ « 2\Ti2s\. The class of NP models 
which do not affect Ti2s then satisfy AF^ < AF™ [11 . 
These include the minimal flavour violating models [l3| 
where the bases in the quark flavor space are the same as 
that in the SM, as well as models where the mixing box 
diagram contains only heavy degrees of freedom. The 
predictions of these models for (AFs, Ps ) will then be 
restricted to the sinusoidal band shown in Fig. [T] Note 
that only a small part of this band is within the 68% C.L. 
region, so that NP of this type will be unable to account 
for the measurements if the errors decrease with the best 
fit values staying unchanged. 

However, there are well-motivated models where the 
Bs-Bs mixing box diagram contains two light degrees of 
freedom, resulting in an absorptive amplitude. Given the 
current strong constraints on the Bs decays to hadrons, 
e~^e~ and ti~^fJ,~ [^, the only candidate for the interme- 
diate light particle is r. In an earlier publication [14| . 
we had implemented this idea with two examples: (i) 
the model with a scalar leptoquark, and (ii) R-parity vi- 
olating supersymmetry. These models can have flavor 
dependent couplings of a light particle with a heavy new 
particle - in particular, r can couple with the leptoquark 
or squark - and hence can contribute to Ti2s- A signifi- 
cant enhancement of AFg was shown to be possible in the 
former model [ij] . In this paper we shall investigate the 
effect of the leptoquark on the correlation between AF^ 
and /3s , keeping in mind that any such new physics 
will also significantly affect the decay rate Bg -^ t~^t~. 

Recently, the D0 collaboration has claimed evidence 
for an anomalous like-sign dimuon charge asymmetry [15| 



A°i = -0.00957 ± 0.00251 ± 0.00146 



(7) 



CDF has also measured the same quantity using 1.6 fb ^ 
of data and found ^^i = (8.0 ± 9.0 ± 6.8) x IQ-^ ^. 
Combining these two, one gets 



At 



.5 ±2.8) X 10" 



which differs from the SM prediction 

,6(SM) _ „ p,r,p,r,o-|-0. 00005 

^sl ~ ~'J-'JU'J^^-0.00006 



(8) 



(9) 



by about 3a. Such an asymmetry can be used as a probe 
of the fiavor structure of new physics [ij]. It turns out 
that the same new physics that enhances AFg can also 
help in explaining this anomaly. We shall elaborate on 
this in the latter part of this paper. 



II. NEW PHYSICS IN Bs -^ t+t' 

Leptoquarks (LQ) are color-triplet objects that couple 
to quarks and leptons. They occur generically in GUTs 



[18| , coniposite models [19|, and superstring- inspired Eq 
models [20|. Model-independent constraints on their 
properties are available j2ll |. and the prospects of their 
discovery at the LHC have also been studied [22] . 

The direct production limits depend on the LQ model, 
as well as the SM fermions these LQs can couple to. The 
bounds on the second and third generation leptoquarks 
are, respectively. My p > 316, 245 GeV, when they are 
pair produced 0, |23|- A third generation scalar lepto- 
quark decaying only into a 6-quark and a r lepton has 
a mass bound of 210 GeV 



24\ 



We shall conservatively 
take Mlq = 250 GeV in this analysis. However our re- 
sults hold even with much higher Mlq, by appropriately 
scaling the coupling |/ilq| as shall be seen later. 

We shall restrict ourselves to scalar leptoquarks that 
are singlets under the SU(2)l gauge group of the SM. 
This is because vector or most of the SU(2)i nonsinglet 
leptoquarks tend to couple directly to neutrinos, hence 
we expect that their couplings are tightly constrained 
from the neutrino mass and mixing data. This makes 
any significant effect on the Bg-Bg system unlikely. 

The relevant interaction term for a scalar SU(2)l sin- 
glet leptoquark is of the form 



-Clq = K^d^jRCiRSo + h.i 



(10) 



where d^ and cr stand for the right-handed down-type 
quarks and right-handed charged leptons, respectively, 
and I, j are generation indices that run from 1 to 3. The 
couplings Xij can in general be complex, and some of 
them may vanish depending on any fiavor symmetries in- 
volved. We take the LQ couplings in the quark mass ba- 
sis. This is the most economical choice given the fact that 
we do not know the rotation matrix for the right-chiral 
down-type quark fields. One can also have an SU(2)i 
doublet leptoquark, whose interaction is of the form 



-Clq = XijqjL^<J2eiRSi + h.c. 



(11) 



which gives almost identical results. 

When A32 and A33 are nonzero, the interaction in 
eq. (ITU)) generates an effective four-fermion {S + P) ^ 
{S + P) interaction leading to 6 — ?> st^t^ . This will 
contribute to Bg-Bg mixing (with t and Sq flowing in- 
side the box), to the leptonic decay Bg -^ t^t~, and 
to the semileptonic decays B -^ XgT^r^. The relevant 
quantity here is the coupling product 



hhQib 



-^ ST T 



^ A09A, 



32^33 



(12) 



One may get a tight constraint on |/ilq| from Bg — s- 
T^T^ . One expects the lifetimes of Bd and Bg to be the 
same in the Standard Model: tbJtb^ = 1.00 ± 0.01 [1|. 
This is certainly true if we assume spectator dominance: 
the decays which do not have a spectator quark con- 
tribute negligibly in the total decay width. Experimen- 
tally, Tg/Td - 1 = (3.6 ± 1.8)% [13. Thus, the branching 
ratio B = BR(i3g -^ t^t~) can be as large as 6-7%. Con- 
sidering the deviations from the naive spectator model. 



which is expected to be small for the Bd-Bg system, one 
may conservatively put the upper bound for B at 10%. 
The value of B is only 0(10"^) in the SM. This decay 
has not been observed, nor is a direct measurement of an 
upper bound on its branching ratio available. A similar 
estimate of ;B ~ 5% is available in [23. If indeed |/ilq| is 
large enough to cause such a significant enhancement in 
B, it is related to B directly through 



B: 



*LQ 



/. 



Mk 



m^ 



128nMl^ 



Ml 



'1-4 



Ml 



9.5% 



^LQl 



^ /250GeVV 



0.3 ; V Mlq 



Ie 



0.250 GeV 



,(13) 



where /b^ is the Bg decay constant. It can be seen that 
for Mlq — 250 GeV, B ~ 10% can accommodate |/ilq| ~ 
0.3. 

We shall show in the next sections that the values of 
/iLQ I allowed by the above analysis can cause significant 



changes in the values of ATg and /3, 
enhance A^; by a sizeable amount 



J/V0 



and can also 
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,7/^0 



where So{Xr) is another Inami-Lim function, and the 
phase space factor is F{t) = 0.64. The details of the 
calculation may be found in [ij]. 

While the next to leading order QCD corrections and 
the l/rub corrections do not affect Mf^ significantly, 
they modify rf^ by - 30% from its LO value ^. The 
QCD corrections are expected to be different for SM and 
LQ operators since the mediating heavy particle for the 
latter case is a color triplet. The l/rrib corrections are 
also expected to differ since the light degrees of freedom 
that flow inside the mixing box are different too. While 
it is desirable to have an idea of these corrections, since 
we are only showing typical results from allowed lepto- 
quark parameters, such corrections can be absorbed by 
just changing the value of Mlq and the phase of /ilq. 
Therefore in our numerical analysis, we use the SM pre- 
dictions for r^l^ d that include the NLO QCD and 1/mb 
corrections, however for T^^s '^^ only use the leading or- 
der contribution. For the sake of clarity, while calculating 
the combined SM and LQ contribution to ri2s, we use 
only the central value of the SM prediction. Including 
the 30% error in the SM prediction will widen the bands 
for our results shown in Fig. [2l 

In the presence of leptoquarks, eqs. (|5I14I15P lead us 
to write the width difference as 



for 



In the presence of NP contribution, the expressions 
the dispersive and absorptive parts of Bg-Bg mixing can 
be written as 



AF., 



2\Tlfg\Rrcos{ 
Ar™i?rcos(( 



W-0r-2/?n 



5A/ 



(20) 



My 



M 



SM 



Ml?g 



M^'' Rmc' 



i 12s - i 12s 



r 



LQ _ T^SM 



Rre'' 



(14) 
(15) 



The standard model contributions, to leading order (LO) 
in l/rub and ag{mb), are given by |26l[27| 



-pSM _ 
^ 12s — 



iVtbV;gf^XB^VB^MlSo{xt) , (16) 

-[(K6v;*j2r- + (K6iCs)'r™ 

+ 2(K6K:KfcK*s)rT , (17) 



where xb^ = ^b, Bb, Jb ' ^^^ ^^'' ' ^^'^ absorptive parts 
of the box diagrams (without the CKM factors) with 
quarks i and j flowing inside the loop, are given in [5J. 
The short distance behavior is contained in 773 , which 
incorporates the QCD corrections, and in the Inami-Lim 
function So{xt)- The value of Ff|^ has been calculated 
up to 0{l/ml) in [2^, wherein some NP contributions 
to ATg have also been studied. 

The leading order leptoquark contributions to the 
above quantities are [ij| 



M- 



LQ 
12s 



"LQ 



384^2^2^ 



.LQ(0) 
12s 



"•LQ 



2567rM^Q 



XBsflB,So{Xr) , 



XBmlF{T) , 



(18) 
(19) 



where the approximation uses j3f^ « 0. The allowed 
values of /ilq permit Rrcos{4>M — 0r) > 1, so that the 
value of AFg can be enhanced in this model. Fig.[2]shows 
that the enhancement can be even up to ATg w 0.4 ps~^ 
for \hLQ\~ 0.3. 

The decay Bg -^ J/tp<p exhibits CP violation through 
the interference of mixing and decay. The CP violat- 
ing phase measured through the time dependent angular 
distribution of this decay is 



^^ ^ 2^'H M- 



12s 



^ p.}/4„P{SM) 



'^,(21) 



where we have used the approximation |Fi2s| <IC |Afi2s|. 
Clearly at low values of I^ilqI, the allowed range of (pM 
will be restricted to be near zero, and hence f3g will 
be close to its SM value, which itself is close to zero. For 
higher |/ilq|, however, the value of (J)m can be anything, 



and hence 



,j/^0 



can be anywhere in its conventional 



range [— 7r/2,7r/2]. This is illustrated in Fig. [51 

Figure [2] overlays our predictions with the leptoquark 
model in the ATs-f3s ' plane on the results of the com- 
bined analysis of CDF and D0. Clearly, the additional 
leptoquark contribution not only can enhance AF^ and 
/3s, but also can allow us to be well within the 68% C.L. 
region of the current best flt. 




JT D0:6.1fb~UcDF: 1.6fb "^ 
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FIG. 2: The predictions of (AF^,^/^''"*) within the scalar 
leptoquark model, overlayed on the combined experimental 
constraints by CDF and D0 through Bs ->■ J/ip4> (Fig- [III- 
Magenta (dark gray), black and aqua (light gray) bands cor- 
respond to |/ilq| — 0.07, 0.17 and 0.27, respectively, with 
Mlq = 250 GeV. 



FIG. 3; The predictions for (Ars,/33') within the scalar lep- 
toquark model, overlayed on the 68% (continuous blue) and 
95% (dashed red) C.L. contours for the combined D0 and 
CDF measurements of afi. Magenta (dark gray), black and 
aqua (light gray) bands correspond to |/ilq| = 0.07, 0.17 and 
0.27, respectively, with Mlq = 250 GeV. 



IV. NEW PHYSICS IN A 



si 



The like-sign dimuon charge asymmetry A^j measured 
by D0 [l^ and CDF [1^| is related to the semileptonic 
decay asymmetries a^j and afj in the Bd and Bg sectors, 
respectively, through [15;] 



All = (0-506 ± 0.043) afi + (0.494 ± 0.043) a; 



si - 



(22) 



The coefficients here are valid even in the presence of NP. 
The average A^j from eq. ([5]), and the current experimen- 
tal constraints of afj = -0.0047 ± 0.0046 US, yield 

<i = -0.012 ±0.007, (23) 

which is almost 2a away from the SM prediction [5| 



s(SM) 



= (2.1 ±0.6) X 10" 



(24) 



This quantity is directly related to AFg and the Bg-Bg 
mixing phase via 



AF, 
AM,, 



■ tan( 



AF, 
AM. 



tan 2/3; 



■si 



(25) 



where (j)f = Arg(— Afi2s/ri2s) — (j)s and we have defined 
I3f such that (j}f — —2(3f. From eq. ([6]), we have 



/3^' = -0.0020 ± 0.0003 



(26) 



In the presence of NP that affects Fi2s, eqs. P^ and 
TS)) yield the relation 



(3f = ^Arg 



12s 



^ ^sl(SM) 



'^M 



r±L + ZL. (27) 



Since /3//'^'^(™) « « /?f ™\ eqs. m and m dearly 



show that (if is in general different from f3. 



J/i"t> 



Note 



that when NP does not affect Fi2s, the value of (f>r van- 
ishes and only then can one say (3f w /3s ■ There- 
fore, it is not recommended to superimpose the parame- 
ter spaces of (AF,,;9^i) and (AF,,/3//'^'^). 

In Fig. [21 we show the constraints in the (AFs,/3^') 
parameter space coming from the ^^j (consequently, a|[) 
measurement in [l^, and (ATg^Pf) predictions at some 
allowed |/ilq| values. It shows that the leptoquark con- 
tribution can give rise to a|j values well within the 95% 
C.L. region of the experimental data. Note that the pre- 
dictions shown in Figs. [21 and [31 correspond to the same 
set of NP parameters. This again illustrates the need to 
clearly differentiate between /3f and (3, 



J/i,4> 



V. SUMMARY AND CONCLUSIONS 

The model with a scalar leptoquark, presented in this 
paper, belongs to the special class of NP models that af- 
fect the absorptive part Fi2s of Bg-Bg mixing. It can 
therefore evade the relation AF^ < APf *^ and can give 
enhanced values of the lifetime difference in Bg-Bg sys- 
tem. The enhancement in AF^ also corresponds to an 
enhancement in the branching ratio BR{Bg -^ t^t^). 

Recent measurements of AFs and (3g by the CDF 
and D0 collaborations exclude the SM prediction to 90% 
C.L.. We illustrate with the example of the scalar lep- 
toquark model that ATg as large as 0.4 ps^^ may be 
achieved, and values in the (AF^, f3s ) parameter space 
close to the best fit from these measurements can be ob- 
tained. Indeed, if future experiments decrease the er- 
rors on these quantities while keeping the best fit values 
at their current positions, only models belonging to this 
class will be able to explain the deviation from the SM. 

The explanation of anomalous like-sign dimuon charge 



asymmetry recently observed at D0 is also facilitated 
by this class of models, since these models give rise to 
large AF^ as well as large (3f simultaneously. We point 
out that these models in general imply that Pf ^ /3s , 
so one has to be careful when including NP in the anal- 
ysis. Also, note that this mechanism affects A^j through 
the modification of AF^ and (f>s , without the need of an 
explicit b — > sfi^ ii~ coupling. This is a common feature 
of all models which have an absorptive part in the Bg-Bg 
mixing diagram. 

In order to confirm the compatibility of such models 
with the data, one needs further NLO calculations of the 
predictions of these models, as well as a better measure- 
ment of Bg — !■ T^T^ branching ratio, which will be cru- 
cial to constrain the leptoquark couplings. The r from 
Bg — > T'^T~ may be expected to have enough energy 
boost at the LHC to be detected. The r polarization can 
also be measured: the r's coming from leptoquarks are 
expected to be right handed. In addition, if we have an 



SU(2)i doublet leptoquark Si , this will also give rise to 
the FCNC top decay t -^ ct^t~ at the level of 1%, which 
will be another probe of the new physics of this class. 
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